Spontaneous avalanche ionization of a strongly blockaded Rydberg gas 
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We report the sudden and spontaneous evolution of an initially correlated gas of repulsively 
interacting Rydberg atoms to an ultracold plasma. Under continuous laser coupling we create a 
Rydberg ensemble in the strong blockade regime, which at longer times undergoes an ionization 
avalanche. By combining optical imaging and ion detection, we access the full information on the 
dynamical evolution of the system, including the rapid increase in the number of ions and a sudden 
depletion of the Rydberg and ground state densities. Rydberg-Rydberg interactions are observed to 
strongly affect the dynamics of plasma formation. Using a coupled rate-equation model to describe 
our data, we extract the average energy of electrons trapped in the plasma, and an effective cross- 
section for ionizing collisions between Rydberg atoms and atoms in low-lying states. Our results 
suggest that the initial correlations of the Rydberg ensemble should persist through the avalanche. 
This would provide the means to overcome disorder-induced-heating, and offer a route to enter new 
strongly-coupled regimes. 



Ultracold plasmas (UCPs) formed by photo-ionizing 
ultracold neutral atomic or molecular gases [l], Q offer 
an ideal laboratory setting to better understand exotic 
phases of matter such as dense astrophysical plasmas Q 
and laser induced plasmas 0|. Experimental and theo- 
retical progress is driven by the potential to reach the 
so-called strongly-coupled regime [1, 0], in which the 
Coulomb interaction energy dominates over the kinetic 
energy of the ions, giving rise to collective effects and 
strong spatial correlations between particles. It is quan- 
tified by the coupling parameter Y = q 2 /4weoakBT 3> 1, 
where a is the Wigner-Seitz radius, q the electron charge, 
and T the ion temperature. In laser-cooled gases, Y ~ 
0.1 — 2 is readily achieved, which has allowed the ob- 
servation and driving of collective mechanical modes 0] ■ 
Reaching deep into the strongly-coupled regime, however, 
has remained out of reach partly due to disorder-induced- 
heating (DIH), where the Coulomb interaction energy 
due to the initially random distribution of the atoms is 
converted into kinetic energy of the ions [H, Q . 



Gases of atoms in high-lying (Rydberg) states offer an 
alternative approach to study UCPs, as they can be eas- 
ily ionized, and the strong Rydberg-Rydberg interactions 
lead to dramatic new effects. The spontaneous evolu- 
tion of an attractively interacting Rydberg gas into an 
UCP has been observed in several experiments 13, U|, 
which initiated in-depth studies on the ionization mecha- 
nisms [l2-[l4| and on the electron- ion recombination dy- 
namics towards Rydberg states 15j. Recently 16j, the 
long-term formation of an ionic cloud from attractive 
and repulsive Rydberg states has been observed, and 
its back-action onto Rydberg laser-excitation rates has 
been characterized. At the high densities achievable 
using optical or magnetic traps, the Rydberg blockade 
effect [1 7| g ives rise to strong correlations in the ini- 
tial gas [LSI]. Since these correlations resemble those in 
strongly coupled plasmas [l9j], Rydberg blockade should 
help mitigate DIH and provide a path to new strongly 
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FIG. 1. (Color online) Generation of an ultracold plasma 
from a Rydberg-blockaded gas. (a) An initially uncorre- 
cted gas of ultracold atoms is prepared in an optical dipole 
trap (I). By coherent laser coupling we excite some atoms to 
a Rydberg state. The strong Rydberg-Rydberg interactions 
prevent the excitation of close-by pairs, leading to spatial cor- 
relations (II). After a time t cf it the Rydberg gas is observed to 
spontaneously ionize into an ultracold plasma (III), (b) Qual- 
itative dynamics of the involved populations: ground state 
atoms N g (solid line), Rydberg atoms N r (dashed line), and 
ions Ni (dash-dotted line). 



coupled regimes |8l. 1 2 Of . However, mechanical collapse 
driven by attractive forces HI, 21 1 and anisotropic inter- 
actions could destroy the correlations. 



In this letter, we report the dynamical formation of an 
ultracold plasma from an initially spatially-correlated gas 
of isotropic and repulsively interacting Rydberg atoms 
in the blockade regime. Despite the relative stability of 
repulsively interacting Rydberg ensembles, we observe, 
after a well defined but controllable time, a spontaneous 
avalanche ionization, as evidenced by a sudden increase 
in the number of detected ions, accompanied by a deple- 
tion of the ground and Rydberg state populations. We 



2 



point out the relevant processes leading to plasma for- 
mation in this system, and using a simple coupled rate- 
equation model we quantify the relevant formation rates. 
We observe that repulsive interactions delay the onset 
of avalanche ionization, leading to a new density depen- 
dence. We estimate that the typical time-scale for plasma 
formation is short compared to the motional dynamics of 
the ions, suggesting that the initial correlations should be 
preserved in the plasma phase. 

The basic principle behind our observations can be un- 
derstood in three main stages (Fig. [T]). (I) An initially 
randomly distributed gas of ultracold neutral atoms is 
excited to Rydberg states via a continuous two-photon 
laser coupling. (II) Due to the Rydberg blockade, each 
Rydberg-excited atom blocks further excitations within 
a radius R c leading to density-density correlations which 
resemble those of a gas of hard-spheres (Figj2]- inset). Af- 
ter a short time, the Rydberg density reaches steady 
state, however over time Rydberg atoms start to decay 
by a combination of blackbody photoionization [l2| and 
ionizing collisions with atoms in the ground and interme- 
diate states 22j , which leads to a gradual increase in the 
number of charged particles in the system. (Ill) Once 
a critical number of ions N CT ix accumulates the resulting 
space charge can trap subsequently produced electrons. 
At this critical time t cr it, rapid electron- Rydberg colli- 
sions trigger an ionization avalanche, leading to the for- 
mation of an UCP. Since the avalanche is triggered by the 
fast moving electrons the plasma forms rapidly compared 
to the motion of the ions. Therefore the original spatially 
ordered structure of the Rydberg ensemble should be pre- 
served. This is predicted to reduce the effects of DIH by 
several orders of magnitude Q. Due to the continuous 
laser coupling the avalanche continues until the ground 
state population of the gas is fully depleted. 

In our experiments we measure simultaneously the 
atom number in the ground or intermediate state N g , 
the Rydberg atom number N r , and the ion number Ni 
in a small excitation volume through a combination of 
absorption imaging and ionization detection. Our exper- 
imental setup is as follows: dense atomic samples of up 
to 1.5 x 10 5 87 Rb atoms in the \5S 1/2 ,F = 2) state are 
loaded into a single beam optical dipole trap. We hold 
the atoms in the trap for a short time to equilibrate at a 
temperature of T = 230 ± 30 fiK. After a typical time-of- 
flight of 200 fis the cloud has expanded to a radially sym- 
metric Gaussian distribution with e _1 / 2 -radii of er ra diai = 
28 ± l^m and cri ong = 310 ± 10/im. The peak atomic den- 
sity is varied between uq = (1 — 10) x 10 10 cm~ 3 while 
maintaining an approximately constant cloud shape by 
reducing the overall atom number during trap loading. 

Rydberg atoms in the |55 i S , 1 / 2 ) state are excited us- 
ing a doubly resonant narrow-band two-photon excita- 
tion via the intermediate I5P3/2) state. Both lasers are 
aligned perpendicular to the symmetry axis of the cloud. 
The first excitation step is realized with a nearly uni- 
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FIG. 2. (Color online) Density dependence of the Ry- 
dberg number for short times. The atoms in \55S) are 
field-ionized after 3.7/j.s of continuous laser coupling. Strong 
interactions lead to a Rydberg excitation suppression for in- 
creasing ground state densities, in contrast to the interaction- 
free case (dashed line) . The data is in good agreement with a 
hard sphere model (solid line) with a sphere diameter corre- 
sponding to the blockade radius R c ~ 5/im. Inset shows the 
spatial correlation function typical for a Rydberg blockaded 
gas El. 



form beam at 780 nm wavelength with intensity I p — 
2.9mW/cm 2 . The second laser, at 480 nm, is aligned 
counter propagating to the first and is focused onto the 
center of the cloud with a Gaussian waist w = 40 /jm and 
a peak coupling strength S1/2-7T = 5.8 MHz. The excita- 
tion volume V, defined by the intersection of the cloud 
and coupling beam, is approximated by a cylinder along 
the coupling beam, of length L = v27r<T ra dial ~ 70 /im 
and radius R = w/\/2 ss 28 /im. To perform the exper- 
iment we simultaneously pulse on both lasers for times 
t ranging from to 38 jj,s. During this time, V changes 
by less than ~ 20%. To stop the excitation, we switch 
on an electric field which field ionizes the Rydberg atoms 
and subsequently guides all ions to a micro-channel plate 
detector (MCP). 

We first investigate the dependence of the Rydberg 
population on the density uq for a short excitation time 
t = 3.7 /is (Fig. [2]), at which we can rule out the influence 
of spontaneously formed ions. We observe the character- 
istic nonlinear dependence and saturation of N r which is 
a consequence of Rydberg blockade 23]. The observed 
scaling is well described using a classical har d sp here 
model in the steady-state (solid line in Fig. [2]) 24]. In 
this model, each Rydberg-excited atom is assumed to 
produce an exclusion sphere of radius R c rj y/ {2C(X)/VL 2 
which reduces the fraction of ground state atoms avail- 
able for subsequent excitation. From the strength of the 
van-der-Waals interactions, C§/2n s» 50 GHz /im 6 for 
the |55S , i/ 2 ) state, and from the intermediate-state decay 
rate, r/27r w 6.1 MHz, we get R c ~ 5 /im. We assume 
independent excitation A x cf> and de-excitation rates B. 
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Interactions are accounted for through the fraction of 
available ground-state atoms 4> < 1, which decreases for 
increasing Rydberg density and is estimated using the 
Carnahan-Starling expression for hard spheres (see Sup- 
plemental Material and [Hj]). The ratio A/(A+B) = 0.38 
is constrained by the steady state solution to the (single- 
atom) three-level optical Bloch equations for our exper- 
imental parameters. Results of the analytic model were 
verified using Monte-Carlo simulations of the excitation 
process. 

From this model we derive the steady state N r as a 
function of n assuming homogeneous density and cou- 
pling strength in V. As seen in Fig. [2] the hard-sphere 
model is in excellent agreement with the data, with the 
overall detection efficiency ~ 0.4 as only adjustable pa- 
rameter. From comparison to the data we find that inter- 
actions suppress N r by up to a factor of 8. At the highest 
densities the fraction of blockaded atoms 1 — <j> ps 0.9, and 
the packing fraction n = ttN t R^/6V ~ 0.2, which indi- 
cates the presence of strong correlations in the system. 

To study the evolution of the system we probe the 
ground state atomic distribution using resonant absorp- 
tion imaging. This is performed using a second light 
pulse on the lower transition, applied 1 (is after the 
end of the Rydberg excitation. Figure^ shows a se- 
quence of absorption images of the center of the cloud, 
recorded for different excitation times, at a density no — 
2.8 x 10 10 cm~ 3 . For times longer than pa 10 /zs we observe 
the appearance of a hole coinciding with the position of 
the focused 480 nm excitation laser. Over 10 — 20 /is 
the hole grows until the cloud becomes locally trans- 
parent. We estimate the remaining number of atoms 
N g in V (shown in Fig. [5b) by fitting to the images the 
difference of two peak functions. The time-dependence 
shows initially a slow decrease of N g , until a critical time 
£ crit rj 12 its at which point starts a rapid depletion. 

This behaviour is a signature of the ionization 
avalanche of the Rydberg population. To independently 
measure N r and Ni, we perform two experiments: in the 
first, we use an electric field above the Rydberg ioniza- 
tion threshold, and we detect both Rydberg atoms and 
ions using the MCP. In a second experiment, we switch 
a field below the ionization threshold, which has a neg- 
ligible effect on the Rydberg atoms but is sufficient to 
guide the ions to the detector. After relative calibra- 
tion of the detection efficiencies for both experiments, 
we substract the two signals to deduce N r and Ni indi- 
vidually (Fig. [3Jc) . N r reaches a quasi-steady-state value 
s» 700 after ps 4/zs. Around t crit ~ 12 us, the number 
of ions reaches JVj = 120 and starts to increase rapidly. 
This number is already sufficient to trap electrons, thus 
leading to the formation of an UCP [6] , with lifetime 
ps 9 /xs. Applying a small electric field during excitation 
suppresses the avalanche ionization, thereby confirming 
that the plasma formation is related to an accumulation 
of charges. 
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FIG. 3. (Color online) Measurement of the complete 
avalanche ionization dynamics for N g , N r and Ni. (a) 

Absorption images taken after different durations of con- 
tinuous laser coupling on a cloud of initial density no = 
2.8 x 10 10 cm -3 . (b) Number of ground-state atoms N g in the 
excitation volume. The data are compared to a rate equa- 
tion model (see text), accounting for (solid line) or neglect- 
ing (dashed line) the Rydberg blockade mechanism. Vertical 
lines highlight the avalanche onset time, (c) Simultaneously 
recorded Rydberg atom population iV r (circles) and ion pop- 
ulation Ni (triangles). 



The observed dynamics can be described by using 
a simple coupled-rate equation model which incorpo- 
rates the Rydberg blockade effect, seed ionization pro- 
cesses, and an avalanche term driven by ionizing electron- 
Rydberg collisions (see Supplemental Material). We ne- 
glect the role of Rydberg-ion interactions, assuming the 
net charge imbalance is relatively small and the ions are 
effectively screened by the mobile electrons. This as- 
sumption is consistent with our experimental observation 
that the Rydberg resonance does not significantly shift 
after plasma formation, but only becomes broadened [l6| . 
On the observed timescale, the repulsive |555' 1 / 2 } atoms 
are not expected to undergo ionizing Rydberg-Rydberg 
collisions at a sufficient rate to seed the avalanche [2o| . 
Instead, we attribute the seed pro cesses to a combination 
of blackbody photoionization 12| and ionizing collisions 
either with ground-state atoms or with the relati vely 
large population in the intermediate state I5P3/2) [22J. 
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Once Ni exceeds iV cr it , we consider N e — N — iV cr it elec- 
trons to be trapped in a plasma, which subsequently in- 
teract with the remaining Rydberg population. 

We fit the model results simultaneously to a set of 
ground state depletion curves, with initial atomic den- 
sities ranging from (0.8 — 9) x 10 10 cm -3 . The black-body 
photoionization rate is taken from ref. [ijj]. This leaves 
three free parameters that we extract from the model: 
the overall Rydberg excitation rate A from the ground 
and intermediate states, the effective seed collisional ion- 
ization cross section a co i 27J], and the average energy of 
trapped electrons E e , which constrains both -/V cr ;t and 
the non- linear rate for electron- Rydberg collisions j av . 

Typical predictions of the model are shown in Fig. [3]d- 
c. The sharp evolution of N g is well reproduced, although 
the model does not completely reproduce N r and Ni after 
the avalanche. We attribute this to the plasma expansion 
and three-body recombination dynamics, not included in 
the model. The best fit parameters and uncertainties are: 
A = (8.7 ±0.3) x 10- 2 MHz; cr col = 0.73 ± 0.21 ^m 2 , a 
factor of 8-15 larger than the geometrical cross-section of 
the Rydberg atoms (a geo « 7m* 4 ag ~ 0.064 ^m 2 , where 
ao is the Bohr radius and n* is the effective principal 
quantum number); and E e = (3.0 ± 0.9) x 10~ 3 eV, 
i.e. approximately 60% of the electron binding energy 
in 1555*1/2), similarly to[ll|. Given our geometry, E e 
relates to -/V cr ;t = 26 ± 7 (typically 5 % of the ini- 
tial Rydberg number), and to 7 av ~ a geo ^/E e /m e sa 
(1500 ± 200) MHz • ^m 3 [H. 

To investigate the regimes in which a plasma develops, 
we extracted £ cr ;t from the data as a function of density 
(Fig. HI). We observe a critical dependence, following an 
approximately power law of t CI it oc a with a = 1.0 ±0.1 
on the considered density range. We have also studied the 
role of the excitation volume by performing experiments 
for a shorter time-of-flight of 100 [is, which leads to a 
reduced volume deduced ~ V/2. We observe that the 
onset times are delayed, since N CT it decreases slower with 
V than the total atom number Q . This effect is well 
reproduced by the model, when accounting for the new 
geometry with conserved rates and cross sections. 

From the model we conclude that the Rydberg block- 
ade has a significant effect on the formation of the UCP, 
delaying f crit by typically 80% for the data shown in 
Fig. |U In such a time, Rydberg interactions are sus- 
pected to develop strong spatial correlations fl8| . Our 
model predicts a sudden avalanche ionization of the com- 
plete Rydberg ensemble within w 0.2ns. Within this 
time, the ions only move by at most 60 nm, which is much 
less than their average spacing (> R c ): thus, the correla- 
tions from the Rydberg sample are most-likely preserved 
in the avalanche. The observed packingfraction n s» 0.2 
is similar to the highest considered in [8j: provided the 
laser coupling are turned off immediately after t cr i t , we 
expect that the importance of DIH as compared to final 
ion interactions should be strongly reduced. 
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FIG. 4. (Color online) Density dependence of the ion- 
ization avalanche onset time tcrit (log-log scale). Dark 
and light symbols correspond to experiments with excitation 
volumes V and Koduced ~ V/2, respectively. Circles and tri- 
angles correspond to the |55£?i/2) and |40Si/2) states, respec- 
tively. Continuous lines are the model predictions for |55Si/2) 
for both volumes, using the same parameters fitted on exper- 
iments in V. The dashed line shows the model outcome in V 
when neglecting the blockade effect. 



To further study the effect of Rydberg interactions 
on the avalanche dynamics, we compared the dynam- 
ics from different Rydberg states, while keeping all other 
experimental parameters (in particular f2) approximately 
constant. In comparison to |55Si/2), the |40Si/2) state 
appears more stable, exhibiting plasma formation in 
^■educed only above densities of no s» 1 x 10 11 (trian- 
gles in Fig. 2]). Rescaling the parameters <J co i,Jav oc n* 4 - 
R c oc n* 11 / 6 ; and E e 



oc 1/n 



*2 



ll| . the model accounts 
for approximately half of the relative shift in t cr it. For 
the I55D3/2) state we observe an immediate production 
of ions (tcrit 0), presumably due to the anisotropy of 
D-state interactions which breaks the blockade effect and 
leads to an increased rate of ionizing Rydberg-Rydberg 
collisions [l3j . 

In conclusion, we have experimentally investigated the 
spontaneous evolution from a strongly blockaded and 
spatially correlated Rydberg gas to an UCP. Simultane- 
ous absorption imaging and ion detection allow us to ob- 
serve the complete dynamics of this process. In particular 
we identified a regime in which the ionization avalanche is 
faster than the motional timescale of the produced ions, 
suggesting a new route to the production of strongly cou- 
pled UCPs with large ionic coupling parameters T. In 
addition, Rydberg states allow for a large tuning of the 
correlation length given by R c oc n* 11 / 6 . Future experi- 
ments could directly resolve the spatial structure of the 
resultant plasmas, for example by using Rydberg states 
as a sensitive probe [2^]. This would open promising 
new avenues to experimentally investigate novel effects 
in UCPs, such as liquid-vapor phase coexistences, criti- 
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cal points [30( and strongly coupled liquid phases. 

This work is supported in part by the Hei- 
delberg Center for Quantum Dynamics and the 
Deutsche Forschungsgemeinschaft under WE2661/10.2; 
M.R.D.S.V. (grant number FP7-PEOPLE-2011-IEF- 
300870) and S.W. (grant number PERG08-GA-2010- 
277017) acknowledge support from the EU Marie-Curie 
program. 



SUPPLEMENTAL MATERIAL 

Coupled rate-equation model of the avalanche 
dynamics 

We describe the avalanche ionization by a coupled rate- 
equation model between the number of atoms in the 
ground or intermediate state, both included in N g , the 
number of Rydberg atoms N r , the number of ions N, and 
the number of electrons N e trapped in the plasma. Only 
two coupled differential equations describe the evolution 
of the system: 



dN, 

dt 
dN r 

-dT = A(pI "° 

-lcoiN r N g /V - lav N r N e /V . 



= —A<f)N g + BN r 

= A<t>N g - BN r - j bbi N r 



(1) 
(2) 



Here, V is the excitation volume and ibM ~ 150/s is the 
black-body photo-ionization rate • 

We deduce N and N e from the conservation of the 
total particle number N tot , neglecting expansion of the 
plasma. Considering the electron trapping criterion 
we get: 



N = Nt. 



N n -N r 



'tot J 'g 

N e = Ni-N ait , if > . 



(3) 
(4) 



The fraction <p of ground state atoms that, despite the 
blockade mechanism, is available for Rydberg excitation, 
is evaluated following the Carnahan-Starling approxima- 
tion for hard-spheres 25| : 



exp 



-877 + 9rj 2 - 

(i-v) 3 



(5) 



where r\ — nN r R^/SV is the packing fraction of hard- 
core spheres with radius i? c /2. We relate the non-linear 
seed ionization rate j co i to an effective collision cross sec- 
tion a C oi between Rydberg atoms and atoms in either the 
ground or the intermediate state: 



led = &col 



16k B T 
TTm Rb 



(6) 



where T is the temperature of the atomic cloud, fcs is 
the Boltzmann constant, and mm, is the mass of 87 Rb. 



The non-linear avalanche rate i av is related to the elec- 
tron energy E e by 



la 



' geo 



(7) 



Here a geo = im^a^ is the geometric cross-section of Ry- 
dberg atoms [28l | , with n* the effective quantum number 
and a,Q the Bohr radius. The electron mass is given by 
m e and we assume that half of E e is of kinetic form. 

Finally, assuming V to be a homogeneous cylinder of 
length L and radius R, we evaluate the critical ion num- 
ber for electrostatic trapping of electrons with total en- 
ergy E e : 



Nr: 



8E e LirR 2 e 



q 2 [L(-L + y/L 2 + 4i? 2 ) + 4i? 2 csch" 1 (2R/L)] ' 

(8) 

where eo is the vacuum dielectric permittivity, q is the 



electron charge, and csch 
cant function. 



is the inverse hyperbolic cose- 
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